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Abstract 

We prove the propagation of regularity, uniformly in time, for the scaled solutions of the 
one- dimensional dissipative Maxwell models introduced in |3J. This result together with the 
weak convergence towards the stationary state proven in [24] implies the strong convergence 
in Sobolev norms and in the L 1 norm towards it depending on the regularity of the initial 
data. As a consequence, the original non scaled solutions are also proved to be convergent 
in L 1 towards the corresponding self-similar homogenous cooling state. The proof is based 
on the (uniform in time) control of the tails of the Fourier transform of the solution, and it 
holds for a large range of values of the mixing parameters. In particular, in the case of the 
one- dimensional inelastic Boltzmann equation, the result does not depend of the degree of 
inelasticity. This generalizes a recent result of Carlen, Carrillo and Carvalho [TT], in which, 
for weak inelasticity, propagation of regularity for the scaled inelastic Boltzmann equation 
was found by means of a precise control of the growth of the Fisher information. 

1 Introduction 

In 2003 Ben-Avraham and coworkers [3] introduced a one-dimensional model of the Boltzmann 
equation, in which binary collision processes are given by arbitrary linear collision rules 

v* — pv + qw, w* — qv + pw\ p > q > 0. (1) 

The positive constants p and q represent the mixing parameters, namely the portion of the 
pre-collisional velocities (v,w) which generate the post-collisional ones (v*,w*). Under the 
hypothesis of constant collision frequency, this mechanism of collision leads to the integro- 
differential equation of Boltzmann type, 

8tf(v,t)= J (jf(v.,t)f(v>.,t) - f{v,t)f(w,t)j dw (2) 

where now (v*,w*) are the pre-collisional velocities that generate the couple (v,w) after the 
interaction and J = p 2 — q 2 is the Jacobian of the transformation of (v,w) into (v*,w*). As 
observed in [3], while in the long-time limit velocity distributions are generically self-similar, 
there is a wide spectrum of possible behaviors. The velocity distributions are characterized by 
algebraic or stretched exponential tails and the corresponding exponents depend sensitively 
on the collision parameters. Interestingly, when there is energy or momentum conservation, 
the behavior is universal. 
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Since the integrals J m u" f(v, t) dv, n > 0, obey a closed hierarchy of equations [4], moments 
can be evaluated recursively, starting from mass conservation. In particular, choosing as initial 
density a normalized probability density fo satisfying 



fo > 0, / f (v) dv = 1 , / vf (v) dv = , fv 2 f (v) 
J Jm Jr 



dv = 1, (3) 



it follows that both mass and momentum are preserved in time, while the second moment 
varies according to the law 

E(t) = f v 2 f(v, t) dv = exp { (p 2 + q 2 - l)t} . (4) 

■Ik 

Special cases include the elastic model (p 2 + q 2 = 1), which is the analogous of the well-known 
Kac model [181 [23] , the inelastic collisions (p + q = 1 ) [I] , the granules model (p + q < 1 ) 
[26], the inelastic Lorenz gas (q = 0,p < 1) [21j . and in addition energy producing models 
(p>l) [30]. 

By @ it follows that the second moment of the solution is not conserved, unless the 
collision parameters satisfy 

2.2 -i 

p +5 =1. 

If this is not the case, the energy can grow to infinity or decrease to zero, depending on 
the sign of p 2 + q 2 — 1. In both cases, however, stationary solutions of finite energy do not 
exist, and the large-time behavior of the system can at best be described by self-similarity 
properties. The standard way to look for self-similarity is to scale the solution according to 
the rule 

t) = y/Elfif (yy/E(¥),t) . 

This scaling implies that J M v 2 g(v, t) dv = 1 for all t > 0. 

The large-time behavior of the density f(v, t) has been studied in [24| . by resorting to the 
Fourier transform version of the Boltzmann equation ([2]), which reads 

ft/tt, t) = fipt *)/(«€. *) - /(£, t). (5) 
which convertes for the scaled solution g(t) into the following 

dtm t) + \ {p 2 + q 2 - 1) t) = m, t) g {pi, t) - *). (6) 

It is worth recognizing from equation ([5]l an equivalent formulation of equation ([J} making 
use of the convolution operator 

&/(«.*) = /**/«(«. *)-/(«.*) 

where we used the shorthand 

'•'•>- J' (i 

The key argument for studying equations © and (|BJ is the use of a metric for probability 
densities with finite and equal moments of order [a], where, as usual, [a] denotes the entire 
part of the real number a: 

^9) = sa P \S^m. (7) 

The metric ([7]) has been introduced in [16] to investigate the trend to equilibrium of the 
solutions to the Boltzmann equation for Maxwell molecules. Further applications of d a can 
be found in [HI [231 IT71 |2H] . 

The study of the time evolution of the d a -metric, with a — 2 + 5, for some suitable 
< 5 < 1, enlightened the range of the mixing parameters for which one can expect that the 
scaled function g(t) converges (weakly) towards a steady profile g^ at an exponential rate. 
More precisely, let us define, for fixed p and q the function 

s P AS) = p 2+s + q 2+s - 1 - ^ (p 2 + g 2 - 1) • (8) 
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Then, the sign of 5 Pl9 determines the asymptotic behavior of the distance d a (gi(t), gi(t)) 
between two scaled solutions gi(t), <?2(t) issued from two initial data /^o, /2,o. In particular, 
it has been proved in [24] that if there exists 8 £ (0, 1) such that S p , q (5) < for < S < S, we 
can conclude that d,2+s{gi (t), <?2(i)) converges exponentially to zero if initially finite according 
to the bound 

<fa+«(fli(t),fla(t)) <exp{-\S p , q (6)\t}d 2 +s(fi,o,f2,o), 5 £ (0,5). 

It has been also proved in [24] that in this case a unique steady state goo exists for equation 
© and for any initial data /o with (2 + S) finite moments, we have for < 8 < S: 

d2+s(g(t),g 00 ) < exp{-\S p ,g{S)\t}d 2 +s(f ,ga O ) -> 0, t^+co. (9) 

On the original non scaled solution f(t), the limit behavior corresponds to a self-similar 

state foo(v,t) = , 1 ffoo I , " I. Note that, by construction, S p . q (0) = 0, and thus 

min^ e ( 0; i){5p, ? } < 0. A numerical evaluation of the region where the minimum of the function 
5p, 9 is negative for p,q £ [0,2] is reported in [24]. This region includes the relevant cases of 
both inelastic (p + q = 1) and elastic (p 2 + q 2 = 1) collisions, as well as the case, among others, 
in which p = q — 1. 

Despite the fact that the large-time behavior of the solution to the Boltzmann-like equation 
([2jl can be described in terms of the d Q -metric, which is equivalent to the weak* convergence 
of measures [13], the strong convergence of the scaled density g(t) towards goo has never been 
proved before. 

A similar problem occurs in dissipative kinetic theory, where the weak convergence of 
the (scaled) solution to the inelastic Boltzmann equation for Maxwell molecules towards 
the homogeneous cooling state with polynomial tails is known to hold [5] in the d Q -metric 
framework, but the strong convergence is still unknown in full generality. A recent paper by 
Carlen, Carrillo and Carvalho [11] shows that in some cases one can prove that the strong 
convergence holds. Their result, however, requires a small inelasticity regime, which in our 
setting of the mixing parameters means p + q — 1, and at the same time 1 — p 2 — q 2 << 1. 
In this case, in fact, one can resort to methods close to the elastic situation, in which the 
(controlled) growth of the Fisher information, coupled with the exponential decay of the 
da-metric allows to prove the uniform propagation of Sobolev regularity, and from this, by 
interpolation, the strong convergence. 

Propagation of Sobolev regularity for both Kac equation and the elastic Boltzmann equa- 
tion for Maxwellian molecules, together with the precise exponential rate of the strong con- 
vergence to the Maxwellian equilibrium M has been proved in [T5] . The advantage of working 
with the classical elastic Boltzmann equation relies on the fact that one can resort to the 
//-theorem. A careful reading of [12], however, allows to conclude that the proof of uniform 
propagation of regularity makes use of the following condition for the distance between the 
solution f(t) and the Maxwellian M: 

Bup|/(&t)-M(0|->O, t^+co. (10) 

While the convergence of H(f)(t) to H(M) implies the L 1 -convergence and therefore (|10[) . 
the same condition continues to hold provided the Fourier transform of the solution to the 
(scaled) Boltzmann equation g(t) and of the stationary state goo satisfy © together with a 
suitable (uniform) decay at infinity in the £ variable, of the type 

(i+Kieiris(£,oi<*. ce« (ii) 

for some positive constants k, K and p. In fact, in this case, for any given R > 0, 

2K 

!$(£,*) - .M£)l < d a (g(t),goo)R a + 
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which implies, optimizing over R, 

m,t)-9oo(0\<C(a,»,K,K)d a (g(t),g oo y / ( a+ » ) , £ e K, t > 0. 

Therefore, in presence of condition (|lip , the decay to zero of the d a -metric implies the decay 
to zero of sup^ gp \g(£,t) — <7oo(£)l for t — > +00. We will remark in Section |H] that condition 
(|11[) on the initial density fo holds provided the square root of fo has some regularity. For 
example, the Fisher information of fo, which controls the ff 1 -norm of the square root, controls 
l£ll/o(OI ( cfr - the Proof in [20]). 

Condition (JTTJ is difficult to prove directly from the equation satisfied by the scaled density 
g(t), due to the presence of the drift term in equation (|6]). To simplify the proof of the various 
bounds we will introduce a semi-implicit discretization of equation (J6j) , that is, for a small 
time interval At, we will consider the solution to 

gfc* + ^)-gfc*) + i (p 2 + g 2 - 1) ^(C, t + At) = g(pi, t)g{& t) - g(£, t). (12) 

Inspired by the integral formulation of the stationary state goo in [8], the solution g(£,t + At) 
to (|12p at time t + At is shown to be a convex combination of the probability densities g(£, t) 
and g(p£,,t)g(q£,,t). Precisely, for p 2 + q 2 < 1 we have 

(At g(T P t,t)g(Tqtt) +(l-At)g(T^,t))-^ TI , (13) 

where 

1 = l-p 2 -g 2 
r 2 

and for p 2 + q 2 > 1 we have 

g(£,t + At) = ~ f (At g(T P Z,t)g(Tqt,t) + (1 - At) g(r^,t)) . (14) 

At J rSf +1 

Eventhough the solution of the semi-implicit discretization is in both dissipative and non 
dissipative case a uniform approximation of the scaled solution g(t) (see Proposition \W\ in 
Section [3]), we have been able to exploit this approximation only in the dissipative case. In 
this case, the integral formulation is also useful to recover regularity properties of the steady 
solution goo. In particular, we shall prove that the steady state has the Gevrey regularity of 
class A, namely 

e" l4|X |fl«,(OI<l, KI>P. 
where the constants A, /j, and p are related both to p, q and to the number of moments which 
are bounded initially. Our main results are summarized into the following statements. In 
what follows we will denote go the initial data of a scaled solution g(t), even if of course 
go — fo, the initial data of the original, non scaled solution f(t). 

Theorem 1 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there is 5 £ (0, 1) for which S P:q (5) < 0, 
for < S < S. Let g(t) be the weak solution of the equation corresponding to the initial 
density go satisfying the normalization conditions and 



\v\ 2+S go(v) dv < +00. 

Jn 

If in addition 



jfi 

JR 



\9°m< (1+mr ki >R > ^ 

for some R > 0, v > and (3 > 0, then there exist p > 0, k > 0, (3' > 0, v' > such that g(t) 
satisfies 

m,t)\<< \ (i6) 
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Theorem [T] is proven in Section [4] The second result is concerned with the regularity of the 
steady state goo and is proven in Section [5] 

Theorem 2 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there exists 5 £ (0, 1) for which 
Sp. q (8) < for < 5 < S so that a non-trivial steady state goo to the Boltzmann equation ((6} 
exists. Let us denote A £ (0, 2) the exponent such that p x + q x — 1 . Then goo is a smooth 
function and belongs to the \-th Gevrey class G A (R), i.e. 

|3oo(0| <e X p(- M |C| A ), |{|>p 
with suitable positive numbers p and p. 

As a byproduct of both proofs of Theorem [T] and Theorem [5] we also get the uniform 
propagation of the Gevrey regularity for solutions g(t) issued from Gevrey initial data. 

Theorem 3 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there exists 5 £ (0, 1) for which 
Sp, q (S) < for < S < S and let us denote A £ (0, 2) the exponent such that p x + q x — 1. Let 
g(t) be the weak solution of the equation corresponding to the initial density go satisfying 
the normalization conditions f5j), and 



I 

J E 



\v\ 2+S go(v) dv < +oo. 



// in addition 

\go(0\<e- m \ \t\>R, 



for some R > 0, v > and j3 > 0, then there exist p > and k > such that g(t) satisfies 

°" <2 , |C|<p, t>o 
" K|e|mm< "' A \ |£|>p, t>o. 



\m,t)\ < 



The previous results are enough to prove the convergence in strong norms towards the steady 
state goo- In consequence of both Theorems [T] and O we can show in fact the uniform in time 
propagation of regularity in Sobolev spaces of high degree 



imi^ ( k)= / i£™ord£ 

with r\ > 0. It is enough to apply the technique developed in [T5] for the Boltzmann equation 
for Maxwell molecules for showing that whenever the equation propagates a tiny degree of 
regularity, as in Theorem[T] this implies that the equation propagates regularity of any degree. 
Then, using the regularity in high Sobolev spaces, we can pass from the weak convergence in 
d a -metric obtained in [53] into convergence in all Sobolev norms, and strong L 1 convergence 
at an explicit exponential rate for a certain class of initial data. This is the objective of 
Section [6] and the main result is summarized as follows. 

Theorem 4 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there exists 6 £ (0, 1) for which 
Sp.q(S) < for < S < 8 and let g^o be the unique stationary solution of ([6]). Let the initial 
density go satisfy the normalization conditions 0), and 

v\ 2+S go(v) dv < +oo. 



// in addition go £ H V (M) for some n > 0, ^/go £ H"(R) for some v > 0, then the solution 
g(t) of converges strongly in L with an exponential rate towards the stationary solution 
goo, i.e., there exist positive constants C and 7 explicitly computable such that 

llflCO -ffoo|Ui(K) < Ce" 7 *, t>0. 
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Thanks to the scaling invariance of the L 1 norm, Theorem [4] allows to deduce also the 
strong convergence of the original non scaled solution f(t) to the self-similar state foo(v,t) = 
1 



ll/(*)-/oo(*)IUi(R) <Ce-*, t>0. 
Finally, in Section [7] we will discuss in details the relevant case in which p + q = 1, which 
corresponds to the one-dimensional inelastic Boltzmann equation for Maxwell molecules pQ. 
In this case, in fact, it is known that an explicit stationary solution to equation © exists, 

3.(1*1) = (l + |£|)e- |fl , 

or, in the physical space 

o 

9oo(v) = 



n{l + v 2 ) 2 ' 

This stationary solution is independent of the values of p and q, and within the set of functions 
/ which satisfy the normalization conditions ©, is the minimum of the convex functional 



H{f) = - f JW)dv. 
Jm 



This suggests the idea that H is an entropy functional for the scaled equation ||B). but the 
proof of this conjecture would require to satisfy an inequality which is the analogous of the 
Shannon entropy power inequality [6] [27], we are not able to prove. 

It is interesting to remark that the ideas of the present paper, which are valid for thep+g = 
1 case, can be fruitfully used for the three-dimensional dissipative Boltzmann equation for 
Maxwell molecules, extending the validity of the recent analysis of [TT] to a general coefficient 
of restitution. However, the proof of the analogous of Theorem [1] to the three-dimensional 
case requires heavy computations, we will publish separately elsewhere. 



2 Preliminary results 

Let us consider the one-dimensional kinetic models of Maxwell-Boltzmann type: 
8 t f(v,t) = J f !/(«„*)/(«;.,*) - f(v,t)f(w,t)) dw 
f(v,Q) = f (v) 



(17) 



where f(v,t) : R x R + — >R denotes the distribution of particles with velocity v € R at the 
time t > and (v*,w*) are the pre-collisional velocities that generate post-collisional (v,w): 

v = pv, + qw-* 
w — qi>t + pWt 

and < q < p. 

Moreover let us suppose fo satisfying the normalization conditions @. In order to avoid 
the presence of the Jacobian we write equation (|17p in weak form, namely 



and 



^ / $(«)/(«, t)dv= II f(v,t)f(w,t) (*(«*)-*(»)) dwdv, 



lim / $(v)f(v,i)dv= / $(v)fo(v)dv 
t— *° Jm Jm 



for any $ bounded and continuous on R. At least formally by choosing $(«) = 1 and 
$(w) = v one shows that, under conditions ((3]) both the mass and momentum are preserved. 
By choosing $(?;) = v 2 we obtain that the energy of the solution 



E(t) = f v 2 f(v,t) dv 
Jm 



G 



satisfies the equality 

E(t) =e (p2+ " 2 - 1)t £(0). 

Therefore, unless p 2 + q 2 — 1, the energy is not preserved. In the dissipative case p + q = 1, 
one has additionally that the moment is always preserved, while the energy is decreasing. 

Following Bobylev [7], the weak form of equation (|17[) is equivalent to the equation in the 
Fourier variables: 

' d t f(£,t) = f(pZ,t)f(qt,t) - f(£,t), t > 

The existence and uniqueness of a solution for any initial data /o satisfying © can be estab- 
lished in the same way as for the elastic Kac equation. 

Theorem 5 (Theorem of existence and uniqueness [241 125) ) 

We consider fo satisfying the normalization conditions £J) and the following Cauchy problem: 
dtf{i, t) = m, t)f(q£, t) - f{t, t), t>0 



Then, there exists a unique nonnegative solution f £ C 1 ([0, +oo), L 1 (R)) to equation f!8\) 
satisfying for all t > 0/ 

/ f(v,t)dv = 1, / f(v,t)v dv = 0. 

JR Ji 

In order to investigate some properties of the solution is useful to introduce the following 
rescaling 



The function g(t) preserves the energy and it satisfies the equation: 

{d t g(Z,t) + ±(p 2 +q 2 -l)td s g(Z,t)=g(pt,t)g(qt,t)-g(Z,t), t>0 

\g(i,o) = MO ==5o(0- 

In the relevant case p + q — 1, equation (|19p admits an explicit stationary state [T] 

5«(0 = (l + l€l)e- |el . 

Pareschi and Toscani proved in [24j that in a quite large range of values of the mixing parame- 
ters, a unique stationary state exists and the unique weak solution g(t) converges to g^. In 
all the cases p + q 7^ 1 (excluding p 2 + q 2 = 1) the stationary state is not explicit, even if some 
properties can be extracted from the analysis of the evolution equation. The convergence 
takes place in a Fourier distance introduced in [16] to investigate the trend to equilibrium of 
the solution to the Boltzmann equation for Maxwellian molecules. We recall the definition of 
this distance. Let < 5 < 1 and let 

M 2 +s = \f>0: I f(v)dv = l, I vf(v)dv = 0, [ v 2 f(v)dv = l, [ \v\ 2+s f(v) dv < 00, 
I Jm Jm Jm Jm 

We introduce on M2+S the distance: 

4w</,p) = ™pJ%iP. 

Let us define, for 8 > 0, 

o /r\ 2+i . 2+5 , 2+6 1 2 . 2 -,\ 

SvA&j^v +q -1 — [p +<? -i). 

The result of Pareschi and Toscani is as follows. 
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Theorem 6 (Pareschi— Toscani [24j) 

Assume < q < p and such that there exists 8 £ (0, 1) for which S Ptq (S) < for < 8 < 8. 
Let g(t) be the weak solution of equation S19\) , corresponding to the initial density go satisfying 
the normalization conditions ^ and 

\v\ 2+S go(v) dv < +oo. 
Then g(t) satisfies for < S < 8 and for c$ > 0: 

\v\ 2+s g{v,t)dv < c s , t>0. 



Moreover, there exists a unique stationary state g x to equation M9[) which satisfies for < 
8 < 8: 

\v\ 2+S goo(v) dv < +oo, 

g{t) converges exponentially fast in Fourier metric towards g^ and the following bound holds 
for0<8< 8: 

d 2 +s(g(t), goo ) < e- |5 ^ W|t d 2+a ( 50 ,3oo). 



3 An iteration process 

The goal of this section is to build up a sequence of functions {g N (£,t)} which approximates 
uniformely the solution g((,,t). In order to do this, for any fixed T > we consider firstly a 
semi-implicit discretization in time of equation (|19|l by partitioning the interval [0, T] into N 
subintervals and we define thus the approximate solution at any time t = j'S for j = 0, . . . , N. 
Secondly, we define g N (£,,t) on the whole interval [0, T] by interpolation and lastly we show 
the convergence of the approximation to the solution. 

In order to lighten the reading, we will postpone almost all the proofs of this technical 
section at the end of the paper in Appendix [8] 



The approximate equation 



In this paragraph we define the approximate solution at any time t = jj^ for j = 0, . . . , N by 
an iteration process and study some of its properties. 

Let T > and At = ^ for TV G N, N > T. Let $f (£), j = 0, . . . , TV be the sequence: 

0o(O = go{0 

0j+l(£)~0f(O 1 > d „jv /*\ . ~N / j-\ a N f >n n at t 

— — — = -Z -j£<Pj+i(0 + Vi mm m) - <p, (£)> j = o,...,n-i 

where - = . 

r 2 

Proposition 7 

Assume < q < p. If go verifies the normalization conditions (J3J), then there exists a unique 
sequence of bounded function tpf for j = 1, . . . , N satisfying \2U\) . 

Proof: Let us begin by proving that 0^ is well defined. In a similar way as in [8] we 
multiply equation (f20|) by (—^7) sgn £ and obtain 

^ (0i (OKI - *) = (-^) sgn £ l^-^- 1 (At 0^{pO0o{qO + (1 - At) ^ it)) ■ 

We assume now p 2 + q 2 < 1. For any positive £ we integrate on [£, +00) and since (po (£) is 
bounded we get: 



0i mr^ = ^ f +a ° (a* 0Z(ps)0Z(q S ) + (1 - At) ^( S )) s-^- 1 



ds. 
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Finally by the change of variables r = s/£ we are led to: 

/■foo . 
(At ^(prO^(qrO + (1 - At) ^(r£)) 



dr 



r + l ' 



(21) 



For any negative £ we integrate on (— oo, £] and in a similar way we obtain that equality (|2ip 
holds for any £ 7^ 0. Moreover since go satisfies the conditions (3) then ipo = go belongs to 
C X (R) and 



ffo(0) = 1 



d 30(C) 



<1 ^f(0) = 0. 



(22) 



Therefore the function ip{ can be defined by continuity in £ = and it is the unique, bounded 
and C X (R) solution of (|2U|) . By an iteration argument the same conclusion holds for any ifif 



obtaining for j = 0, 



£f+i(0 



,N-1 



At 



(At ^(prO^iqrO + (1 - At) <fi? (r£)) 



dr 



For p 2 + q 2 > 1, we repeat the same argument by integrating on [0, f] and [£, 0] and we get in 
the end 



^ AT 



= £ (At 0f (prO#VO + (1 - At) 0f (t£)) 



Applying Fubini's theorem we can remark that for p 2 + q 2 < 1 any ipj' +1 (t;) is the Fourier 
transform of ipj' +1 (v) where for j = 0, . . . , N — 1 



<A) (v) = go(v) 

N 1 \ r 
^+ l(w) = At A 



Ai 



(23) 



with (pf p (v) = (pi an d similarly for yj^. Analogously, for p 2 + g 2 > 1, </3j+i(£) is the 
Fourier transform of ipf + i{v) where for j = 0, . . . , N — 1 



</>o («) = go{v) 

V 

At 



N 1 \ 



[At i * ( 



(1 - At) -tpj 



dr 



(24) 



r+l ' 



In what follows, the function S Pzq (8) is defined as in ©. 
Proposition 8 

Assume < q < p such that there exists S € (0,1) for which S p , q (5) < for < 5 < S. 



Let tfij , for j = 0, . . . , N defined as in 



(|24p . with go satisfying the normalization 



conditions © and J R \v\ 2+s go(v) dv < +00. Then, for < 8 < 6 there exists Cs > such 
that for N large enough (depending on S, T, p and q), for j — 0, . . . , N we get 



<fij (v) > 0, / tpj (v)dv = 1, / (v)dv = 0, / u («) dt> = 1, / |u| yjj (u) dv < C5 

(25) 



Remark 9 

T/ie equalities in imply that there exists C > such that 



mm < i, 



d^f(C) 



d£ 



< C and 



d^f(f) 



d£ 2 



< 1 



(26) 



for any £ £ K , /or iV large enough and for j = 0, . . . ,N. The first two inequalities have 
already been proved in the proof of Proposition^ 
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Definition of the sequence {g N (£,t)} 



N 



We are now in position to define the sequence {g N (£, t)} N - We define g N (£,jji) = ftf (£) for 
j = 0, . . . , N. We extend the definition on the whole interval by interpolation. More precisely, 
let us define: 



where for < t < T we have (Kn — < t < K^jj for Kn € {1, • • • , N} and more precisely 
there is a function < a(t) < 1 such that t = a(t)(K N - l)% + (1 - a(t))K N £. Any g N (£,t) 
is continuous on R x [0, T] and for any t G [0, T] it belongs to C 2 (R). 

The result of convergence is therefore as follows. 
Proposition 10 

There is a subsequence {g N '(£,t)} of {g N (£,t)} which converges uniformely on any compact 
set o/R x [0,T] to the solution g(^,t). 

4 Propagation of regularity 

In this section we prove Theorem [1] Thanks to the uniform convergence of a subsequence of 
the approximate solutions g N {(,,t) to the solution g(£,t) and to the definition of g N (£,t), it 
is enough to prove the bounds (|16p for any ifff (£) uniformly for N € N and j = 0, . . . , N. 
The control of low frequences is a direct consequence of the properties ^ of the initial data 
and of the convergence to zero of the distance d2+6(g(t),goo) and it is proven in the following 
lemma. 

Lemma 11 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there is 8 € (0, 1) for which S p>q (8) < 
for < S < S. Let g(t) be the weak solution of the equation (6$), corresponding to the initial 
density go satisfying the normalization conditions l&j), and 




'K N -1 




t = 




Let <pf the approximation defined in (|20|) . For any < k < 
any fixed T > and any N £ N large enough we get 



| there exists p > such that for 



i 



el < P, t>o, 



(27) 



1 + fc£ 2 ' 



\<Pi (01 < Y^72> \t\^P, 3=0,.-.,N. 



(28) 



Proof: Let us begin by proving (|27[) . We remark first that since for all t > 




1, 




we can deduce at once 



and so for t > and < k < \ there exists p = p(t) > such that 



m,t)\<Y^, \t\<p® 
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but this is not enough because we need an estimate independent of t. We have therefore to 
proceed differently. By Theorem [S] for any fixed < 5 < 8 we have 



d 2+s (g(t), gao ) = sup iMil^MOl < e-I^WI^+aCsoffoo), t > 0. 
Moreover, since > 0, /„ ffoo(«) du = 1, L« ffoo(u) du = 0, J R v 2 goo{v) dv = we have 



so we get uniformly for t > and £ 7^ 0: 

\§i£,t)-§°o(Z)\i,M 

lei 2 



ig(^*)i< + '^'v~r wi Ki ! 



< 1 _ L _,_ (|2) + C , e -|S J ,, g («)|t|^|2+« ) ^ _^ Q 

This shows that for any < k < | there exists p > such that for all t > 

or, which is equivalent, for any < k < 5 there exists p > such that for alH > 

In order to prove (|28[). we would like to exploit again the 0(2+5 distance. For N £ N and 
j = 1, . . . , N — 1 let us estimate first ^2+5(^^+11 <?oo)- We recall that 

0= ^43oo(0 + ffoo(p£)^<»(«£) (29) 



with i = — p 2 9 . By considering the two equations f|20p and (|29p we have 

yf+i(0-g°°(0-(^ y (0-g°°(Q) l./d^ d„ \ 
At = r C ( v df^ +l(5) - d£ 5 °° (5) ) 

In the same way as in Proposition [7] we obtain the integral form: 

0?+iiO - <?°°(0 = ^ J At (ff (P T £)^f (qr£) ~ ff<x,(pr£)&x,(gT£)) 
+ (1-At) (^(rO-SocCrO)-^. 

Therefore for £ 7^ we get: 

i$+i(0-soc(i)i 



|£| 2 + 4 

. _r_ f +ao f A M (prQ^igrQ - goo (prfl&o (grfl | A ^ \0?(tQ - 9oo (tQ\ \ dr 

At A l A ' |e| 2 + a +U ^ |?| 2 +* ir^ + 



Now for t/0 

\0? {pt£)<P? (<?t£) - ffoo(pr0ff<x,(5r|)| 



^ , / JV w 2+b , 2+b\ 2 

< d 2+ sWj ,gao)(p + q )T 



and so 

l^-i(0-5oo(0l ^ r , fN , / /.„ 2+4 , 2+* 



2+<S 



< d2+«(^. 



AT 



5oo ) ^ (Ai(p 2 +* + g 2+a - 1) + 1) ^ 



dr 
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For TV > 112+11 we get 



and remembering that 



S P , q (5) = ( P 2+s + q 2+s - 1) + ^±*(1 - P 2 g 2 ) 



we get 

t (At(p 2+5 + g 2+ * - 1) + 1 

A* - (2 + 5) V 

= At(p 2+S + g 2+a - 1) + 1 
~ 1_ At2±A (l-p 2 -q 2 ) 

_ 1 + S p , q {8) At + (Atf^jl -p 2 - q 2 )(p 2+S + q 2+S - 1) 
l-(At) 2 (2±A(l-p 2 -^)) 2 

Since S p , q (S) < and 

= l+o(At), At-»0, 



, 2 



l_(A^(i±2(l-p2_ 9 2)) 
we get, for TV £ N large enough and j = 1, . . . , TV — 1 



Recursively, we get 



3 H 

where d2+,5(<?o, floo) depends only on the quantities in (|25[1 . Therefore, for j = 1, . . . , TV and 
^/Owe get 

1^(01 < imoi + \t\ 2+6 d2+s(^ , ffoo) < isco(oi + \a\ 2+6 c 

for C = d2+s(go, floo) independent of TV and j. Since 

we get that for any < k < | there exists p > such that for j — 0, . . . , TV 

|pf(e)l<^ \s\<p 

which achieves the proof. □ 

We are now in position to prove Theorem [T] 
Theorem Q] 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there is 8 £ (0, 1) for which S p , q (S) < 0, 
for < S < S. Let g(t) be the weak solution of the equation {jSp, corresponding to the initial 
density go satisfying the normalization conditions and 



\v\ 2+S go(v) dv < +oo. 



If in addition 



\go(0\ < {1 + m y > Kl > R > 

for some R > 0, v > and /3 > 0, i/ien tftere exist p > 0, k > 0, f3' > 0, i/ > s«c/i t/iai g(t) 
|<Kf,*)l< \ \ (® 
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Proof: The bound on the low frequences |£| < p has been established in Lemma 1111 
Moreover, as a consequence of Proposition 2.4 in [T5], we can suppose that condition (|15p 
holds for any |£| > p with a possibly smaller exponent v' . We will prove that for any N £ N 
and j — 0, . . . , N we get 

for positive v' and ft' small enough. By induction we have only to check the bound on 

/ + oo 1 
[At So(prOft(?T£) + (1 - At) g {rt)\ At 

Let |£| > p and r > 1. We are faced with three different situations. 
Case I: If \qr£\ > p, then 

IA.(pr«fc(ffrO| < {l + p pAir 



" (l + /3( P + g)|C|) 



for /3' < /9(p + g). 
Case II: If |prf | > p and |gi~£| < p, then 

\9o(pr^)go(qr^)\ < 

< 



(l+j8prfc|)"' (l + kq 2 T 2 \£\ 2 ) 

1 1 



(i + pp\i\Y' (i + fcg 2 i^i 2 )' 

By choosing v' > small enough we can show that 



1 1 

< 



l 

i +/?!$! y 



Indeed, since ^tg X x < - for any s > we obtain 



< 



i+/3piei; (i + vici 2 ) - \pJ (i + vici 2 ) 

iv' i 



v p/ (1 + kq 2 p 2 ) 

Finally the last term is smaller than 1 for v' < logi(l + kq 2 p 2 ). 

p 

Case III: If \pr(\ < p, then 

\go(pT^)go(qrO\ < 

< 



(l + kp 2 T 2 \£\ 2 ) (1 + kq 2 T 2 \£\ 2 ) 

1 1 



(l + fc P 2 |e| 2 ) (i + fcg 2 |C| 2 ) 

and we have to show that 

1 



< 



(l + fcp 2 |C| 2 )(l + fcg 2 |e| 2 ) " 
Since 

< 



(i + fc p 2|C| 2 )(i + fcg 2 |CI 2 ) " (i + fc(p 2 + g 2 )l?l 2 ) 
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and p 2 + q 2 — C > 0, in order to establish the desired estimate we have only to prove 
that for v' small enough we have 

1 , 1 



or 

i/log(l + / 3|£|)<log(l + fcC|£| 2 ) 
for any |£ > p. This is true since the function 

log(l + fcC|£| 2 ) 



F(0 



log(l + /3|ei 



is bounded below by a positive constant on |£| > p and so by choosing v' > smaller 
than this constant we get the desired estimate. 



Remark 12 

The bounds in ()16p are also equivalent to 

for some positive constants k, p and C suitably chosen. 
Remark 13 

It wouldn't have been possible to prove also the behavior of the solution g(t) on the low fre- 
quencies by induction and this is why we had to exploit the d2+s convergence and the behavior 
of the steady state as in Lerama [TT\ 



5 Smoothness of the steady state 

In this section we prove the following Theorem [2] 
Theorem [2] 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there is 5 £ (0, 1) for which S Pl q(S) < 
for < S < S, so that a non-trivial steady state to the Boltzmann equation ([6| exists. Let 
us denote A G (0, 2) the exponent such that p x + q x = 1. Then g x is a smooth function and 
belongs to the \-th Gevrey class G A (R), i.e. 

<exp(- M |£| A ), \t\>P 
with suitable positive numbers p and p. 

Proof: We are going to prove that g^, G G A (R), following a scheme already used in [22]. 

As in the proof of Lemma \TT\ we know the behavior of g^ close to zero. In particular, by 
the same reasoning, there exist k G (0, i) and p > such that 

IMOI <e~ H \ \i\<P- 

As a consequence of Proposition 2.4 in [T5] and without any loss in generality, we can also 
suppose p > 1. Let us remember that goo satisfies the equation 

+ q 2 - = g^g^qO - fcott) 

and so, denoting again r = 1 ^ a and following [8], we get 

M0 = r/" ^^r (gr ° dr, CGR. 
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Let us consider the space 

X p = {V 6 L°°(R), |V>(£)I < 1. lK0 = M£) for |C| < p} 

endowed with the metric d,2+s defined in ([TJl , where 5 satisfies the assumptions of Theorem [2] 
The space X p is a Frechet space. Let us consider then the function R defined by 

We are going to prove that i? is a contraction on X p and so is its unique fixed point, ft 
is straightforward that R : X p — > X p . As for the contractiveness, for ij) and ip £ X p and 
|£| > p we have 

— |2T5 -l r 7 7^+1 dT(i2 + s WW 



We remark that the assumption 



<[r{p 2+s + q 2+s ) I J^)d 2+ sW,<P)- 



implies 



and so, 



r, /ys 2+5 . 2+« -i,2 + <5 2 2n . r, 

SpA s )=P +Q -H — (1-P -?)<0 



— — (1 -p -q ) < 1- (p ^ + g ^ ) < f 



which is precisely 
So we get 



5 / 1 2 2-, , 2 . 2 



1 - 5 > 1. 



r ( v 2 + s + a 2 + s ) 
d 2+6 (R(iP),R(<p)) < \_2_ s , d2+s(^,<P)- 

2 



Remembering that r = 1 _ p t_ 2 , we get 
and this allows to conclude. Choosing € for example as 



MO 

and defining by induction for n > 



MO iei<p 

|£|>P 



mo lei < p 

we get authomatically d,2+s{ipn, <?oo) — ► for n — > +00, which implies convergence pointwise. 
We show now that there exists /i > such that for all n £ N 

|lM0l<e-" l€|A , \t\>P 

and this uniform estimate passes therefore to the limit and allows to conclude. The only thing 
to control is that for r > 1 and |£| > p we have 

We distinguish three cases: 
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Case I: If qr\£\ > p, since p x + q x — 1 and p > 1 we get 

e^« |A • \MptO\ ■ \MqrO\ < e^^ 1 -^^ < e^^ 1 ^ < 1. 
Case II: If pr|£| < p, then |£| > 1 implies £ 2 > |£] A . Denoting p 2 + q 2 = C, we conclude 

e" l5|A • \MptO\ ■ \MirO\ < e^-^Z 2 ^ 2 ^ < e H«l A -^l«l A (p 2 + 9 2 ) < e l«l A (M- fe C) < 1( 
provided that jj, < kC. 

Case III: Now assume that qr|£| < p while pr\^\ > p. Using the condition p x + q x = 1 once again, 
one finds 

Ml«l A . U/,„C^tM . l./,„C™-tM <- MlSI A -fcr 2 9 2 « 2 -, 1 r A P A |«| A 



e wl • hM<7rf)| • |^(prOI < e 



< e l«| A (M(l-p A )-fc9 2 ) < gl«| A (M9 A -fe^) < 1 



provided that p < 



Remark 14 

_Bj/ i/ie proof of Theorem^ the condition £ G A (R) seems to be sharp, since in the estimate 
of Case I, it wouldn't have been possible to replace A by a > A. Moreover, in the case p + q = 1 
the explicit stationary state is floo(£) = (1 + an d this proves sharpness at least in this 

special case. It is worth noticing that Bobylev and Cercignam in (Theorem 5.3) proved 
that for p + q > 1 and p 2 + q 2 < 1 the stationary state goo satisfies the bounds 

e'^ <l3oo(C)l<(l + !CI)e H51 - 

For the values of p and q for which in addition there is 8 £ (0, 1) such that S p , q (5) < for 
< S < 8, our result improves the upper bound and gives also a new result for some p and q 
in the range p + q < 1 . 

By a careful reading of both proofs of Theorem [T] and Theorem [5] we can deduce that not 
only polynomial tails of the Fourier transform of the initial data go are uniformly propagated 
by the solution g(t), but also exponential ones, as long as the exponent does not exceed the 
exponent A characterizing the mixing parameters p and q. More precisely, the result is as 
follows. 

Theorem [3] 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there is 5 £ (0, 1) for which S Pl q(S) < 
for < 5 < 5. Let us denote A £ (0,2) the exponent such that p x + q x — 1. Let g(t) be 
the weak solution of the equation 0), corresponding to the initial density go satisfying the 
normalization conditions and 

/ \v\ 2+S go(v) dv < +co. 
Jr 

If in addition 

\9o(0\<e- m ", \t\>R, 
for some R > 0, v > and j3 > 0, then there exist p > and k > 0, such that g(t) satisfies 

' V K « 2 , |£|<p, t>0 



\m,t)\ < 



e- K|5|mm "" A) , 151 >P, t>0. 
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Proof: The proof follows the same lines as the proof of Theorem [TJ replacing the polynomial 
decreasing by the exponential one. The key argument for the low frequencies is that the 
normalization assumptions Q on the steady state imply 



which means 



or equivalently 

|&o(0l < e^ 4 *, \i\<p 
for any k £ (0, ^) and p depending on fc. The whole proof follows then without any particular 
difficulty, exploiting for the high frequencies the estimates performed in the proof of the 
Gevrey regularity of the steady state. □ 



6 Strong convergence 

In this section, we are going to prove Theorem [4] on the strong L 1 convergence of the scaled 
solution g(t) to the stationary state goo- 

Theorem [4] 

Assume < q < p satisfying p 2 + q 2 < 1 and such that there exists 8 £ (0, 1) for which 
Sp.q(S) < for < S < 5 and let be the unique stationary solution of ©. Let the initial 
density go satisfy the normalization conditions and 

\v\ 2+S go(v) dv < +oo. 

// in addition go G H V (R) for some n > 0, i/go £ H"(M.) for some v > 0, then the solution 
g(t) of fSJ) converges strongly in L with an exponential rate towards the stationary solution 
goo, i.e., there exist positive constants C and 7 explicitly computable such that 

\\g{t) -SoolUitR) < Ce~ 7t , t>0. 

Let us begin by the following lemma. 

Lemma 15 

Let the initial density go satisfy the normalization conditions and 



\v\ 2+S go(v) dv < +00. 



If in addition ^/go € H"(M.) for some v > 0, then go satisfies 

for positive constants C and (3 and the solution g(t) of fj5j) satisfies 

sup|g(£,i) -300(C)! <C ie - C2t , t>0 (31) 

for positive constants G\ and C2 . 

Proof: Since go = yfgoy/go, then go = yfgo * ^/go- So, for £ £ R we get 

|eriffo(6|< / \i\ v \7^{.i-r)^(j)\ dr 
Jm 1 1 



<K f (|f-T|" + |Tn|v5oK-r)V5o(T)| dr 
7e 1 1 

< 2Jf|| > /floll*" 
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Since moreover |go(f)l — 1; we can nn d positive C and j3 such that 
Thanks to Theorem [1] and to Remark 1121 we get 

\mt)\< J ^ w , ^m, t>o 

for suitable C, n, and \i. The steady state goo belongs to a Gevrey class, so it satisfies an 
analogous estimate, with suitable constants which we can suppose to be the same. Let now 
R > to be chosen in a moment. We get, for (6l: 

2+5 . 2C 



\3{Z,t) -$oo(0l < d 2+ s(g(t), goo )R' + ° + 



{kRY ' 

which implies, optimizing over R, 

|3(^*)-3o (e)|<Cid2 +5 (5(t),ff oo r /(2+5+M) =Cie-° 2t , ^R,t>0. 

for C\ and C2 positive constants. □ 

Proof of Theorem [4] The result consists in converting the weak convergence in the Fourier 
distance of the solution g(t) to the stationary state goo (Theorem [6]) into a L 1 convergence 
by interpolating this weak distance d-2+s with the uniformly boundedness in time of suitable 
moment and Sobolev norm of the solution itself. The only missing ingredient at this point is 
the boundedness of the Sobolev norm and we will go through the proof of it in a moment. 
Let us recall first how we can interpolate these results, following the scheme introduced in 
[12] and fruitfully applied afterward in several papers ([5], [TU] for instance). 

First of all, it is easy to prove the two following interpolation bounds (see Theorems 4.1 
and 4.2 in [E]): for 8 6 (0, 5), there exists a positive constant C such that 

1 2(2+a) 

and for any a > there exist positive constants M, N, f3 and 7 such that 





\Ho\ \ 



l pjpr I (II^IIhm + II^IIh«) 7 . (32) 

So, letting h = g(t) — goa , and s = we get 

- ffoolLi < c (\\\v\ 2+s g(t)\\Li + IIM a+, a»IUi) a x 

d 2+ s(g(t),goof (\\g(t)\\ H M + \\g(t)\\ 

for suitable exponents a, j3, 7. Concerning the stationary state goo, it have been proved in 
Theorem [6] that | v\ goo \\ l 1 ^ for 5 G (0,(5) and thanks to the Gevrey regularity, we 
have goo £ H S (R) for all s > 0. As for the scaled solution g(t), the uniform boundedness of 
the (2 + <5)-th moment have been also proved in Theorem [H] so we will get the L 1 exponential 
convergence as soon as we prove the uniform boundedness in time of g(t) in a suitable Sobolev 
space H max ( M ' N) (R). Of course, we need to assume g e H max( - M ' N) (R). As a byproduct of 
the uniform boundedness of Sobolev norms, we will also get from (|32[1 the convergence of g(t) 
to goo in Sobolev spaces. 

Let us recall how to prove the uniform boundedness of g(t) in a generic homogenous 
Sobolev space H v for r\ > under the assumption 30 £ H v . First of all, let us remark that 
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g(t) £ H v for all t > 0, without any uniformity in time. Indeed, coming back to the original 
non scaled solution f(t) we get 



2 



= 2^ |e| 2 "Re (/(£,*) (/(?£,*)/(«£,*) - /(£,*))) d£ 

= -2 / |e| 2 "l.m,t)! 2 d? + 2 f \e v Re(Wtjf(pt,t)f(qt,t)\ d£ 

<-i/ Ki a "i/«,*)i 2 <«+ / i*i a "i/(pe,«)i 3 i/(tf,f)i a d* 

< -i //H/U *)i 2 de + + ^) jT "leriMt)!- de. 

Since 3 I e 2 '£,*]=/(£, t) we get 
and so 



;llfl(*)lll„ < 



df 

which leads to 



d„ , ll2 

(33) 



llsWII^ = C P , q , v \\g(t)\\%^, 



\\9(t)\\%n < l|<ftl|*,e t W. 



Since p 2 + q 2 < 1, it is not difficult to be convinced that C P , ?I ^ > as soon as for example q 
is small enough. 

Let us make estimate (|33|) more accurate. The goal is to get for example the following 
differential inequality: for two positive constants H and K and to > 0: 

^llffWIll, <-H\\g(t)\\%+K, t>t (34) 

so that 

Let us come back to inequality 
d 



\\g(t)\\%v < Cmax(||g(t )|| 2 j„, 1), t > t . 



ii/wii^<-i ieri/«,*)r<«+ \^r\m,tr\m,t)\^ 



df 



which reads, on the scaled solution g(t), 



f e^^^\\gim 



dt 

<_ e -4^( 2l7+ i )f r ievmt)l2 ^ + e ^^ i 2r l+1)t f \^ lm ,t)\ 2 \g(qtt)\ 2 d^ 



and so 

< (-1- 1 '^~ g2 (2r ? + l)^) ||gW|||,+^|C| 2 "l3(pC,i)| 2 !5(9C,i)! 2 dC. 
Since p 2 + g 2 < 1, it would be enough to obtain for example the following inequality 

\^\ 2v \m,t)\ 2 \9(Q^t)\ 2 ^<^\\g(t)\\ 2 H,+K, t>t (35) 
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where K > is independent of t. Let us prove inequality ((35}. We split the integral in (|35[) 
into two parts 



/ \z\ 2ri \m,t)\ 2 \m,t)\ 2 az= f + f 

Jw. J\e\<R -I i. a 



A + B 

ei>j 



where R will be chosen later. Let us estimate first the term in A (we will denote s a constant 
which is allowed to vary from one line to another, depending at most on p and q). Since 
|<?(£> t)| < 1 for ^ G R and t>0we simply get 

|£| 2 "|3(pU)| 2 |5(9£,;)| 2 d£ 

t\<R 

< I \^^ = -^— i? 2 " +1 , t>0. 

Let us come to the term in B, where we are going to exploit Lemma [15] We remark that 
5oo(£) — > for £ — * +oo and so, by Lemma Il5l for any e > there exist R > and to 
depending on e and p such that 

ifl(pe,*)l<iff(pe,*)-floo(pe)i + i3o (p^)i<2e, |£|>k, t>t . 

We can deduce for t > to: 

z\ 2v \m,t)\ 2 \m,t)\ 2 dt<(2e) 2 [ \e v \m,t)\ 2 ^ 

S\>R • / I€I>-R 

<3^r / i?i 2,7 iff(e^)i 2 dc = eiiffWii^. 

We have obtained 

\Z\ 2v \m,t)\ 2 \g(tf,t)\ 2 dt < £\\g(t)\\ 2 H, + ^p[ R2v+1 ' 1 ^ *°- 

Letting e be fixed such that e < h, we get the desired estimate. 



7 Lyapunov functionals and open questions 

The case p + q — 1 separates in a natural way from the others. It corresponds to a one- 
dimensional dissipative Boltzmann equation in which the momentum is preserved in a micro- 
scopic collision of type ([l}. Equation ([2]) with p + q = 1 as been intensively studied in a series 
of papers [TJ [5] [5] [53] . Among other properties, this model possesses an explicit self similar 
solution, which has been first discovered in pQ. In fact, condition p + q — 1 implies 

\ {p 2 + q 2 - i) = -pq- 

Hence the Fourier transformed version of the scaled equation ([6} can be written as 

9(pZ,t)g(qZ,t)-g(Z,t) = dtg(tt)-pqZd e g(Z,t). (36) 

The choice 

?oo(£) = (l + l?IK kt| 

leads to 

<Mp£)<My'£) - g°o(0 = p<2l£l 2 e~ lel = -p<?£d?i?oo(£) 

and so solves (|36p as a stationary solution for any choice of the parameters p and 5 such 
that p + q = 1. It can be easily verified that in physical variables the steady solution reads 

So» = I . (37) 



20 



Note that this function satisfies the normalization conditions ([3]). Under these constraints, 
however, it can be shown that is the (unique) minimizer of the convex functional 



H(f) = - / ^f(v)dv. (38) 



It is a natural question to investigate whether the functional H is a Lyapunov functional for 
the scaled equation for g(v,t), which can be formally written as 

dtg(v,t) = g p *g q (v,t) - g(v,t) + i (p 2 + q - l) d v (vg(v,t)). (39) 

The results of both Section [4] and Section [6] lead to conclude that, under suitable regularity 
assumptions on the initial value, one can study the time derivative of the functional H(g)(t) 
along solutions to equation (|39|l . obtaining 

d£ - [Jr y/g{v,t) 4 J R 

Hence, the functional H is a Lyapunov functional for equation (|39[) provided the inequality 

i+P 2 + g 2 f /-7-r, ^ f gp*g q {v) 



g(v)dv< I g "\' dv, p + q = l, (40) 

is verified for all functions satisfying constraints ©. We remark that inequality (|40p is 
saturated by choosing g = g x , with g^ defined as in (|37|l . To our knowledge, this inequality 
has never been investigated before, but it can be conjectured that it holds true, even if we 
are not able to prove it. 

A different way to attach the problem is to resort to the Fourier version of equation 
(|39p . This idea has been fruitfully employed in [5] to recover Lyapunov functionals for the 
Boltzmann equation for Maxwell molecules. Let us consider the approximate solution (|13p 
which is a convex combination of the probability densities g(£,t) and g(p^,t)g(q( i ,t). For any 
convex functional H acting on g we obtain 

H(g(t,t+At)) < Jl j °° (At H (g(r P tt)g(Tqtt)) + (1 - At) H(g(rt t))) (41) 

If H{g) = H(g) is defined by (|38]), we have 

H(g(rO) = V^H(g(0), 
and this implies that inequality (|41|l becomes 

Htf{t, t + At))< (At H {g{pi, t)g{qi, t)) + (1 - At) t))) . 

At 2 

Let us suppose that there exists A(p, q) such that for all functions g satisfying conditions (J3J) 
we get 

H (3(PO$(«0) < A(p, q)H(g(0). (42) 

Then, since 



the condition 



would imply 



1 _ l-p 2 -q 2 
r 2 

A(p, 9 ) = i±4±i! (43) 



(AtA(p,q) + (1-At)) = l. 



At 2 

In this case, inequality (|4ip would imply 



H(g^,t + At)) < H(g(i,t)), 
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and H would be a Lyapunov functional. Note that inequality (|42[) corresponds to a reverse 
Young inequality first derived by Leindler [T5]: for < a,/3, p < 1 and /, g non- negative [2] 

ll/*sll/>>ll/ll«llsll/>, 1/a + l/p = 1 + 1/p. (44) 
In our case, p = a = 1/2, f3 = 1 together with the second condition in (|3} implies 

ll/p * /<j||l/2 > P||/||l/2, 

namely inequality (|42p with A(p, q) — p. Unlikely, the direct application of inequality (|44[) is 
not enough to obtain (|42[) . It remains an open question to prove that, under constraints Q 
it holds the Young-type reverse inequality 

\\fp*f4\',l>A{jp,q)\\f\\\',l 
where p + q = 1 and A(p, q) is given by (|43|) . 



8 Appendix 



Proof of Proposition [HJ We will consider only the dissipative case p 2 + q 2 < 1, since the 
other case adapts straightforwardly. It is easy to get for all N £ N and j = 0, . . . , N — 1: 



¥>j+i{v)>0, / tp? +1 (v) dv = l / vipf +1 (v) dv = 0. 



(45) 



Let us consider therefore the evolution of the two other moments of the sequence if j . First 
of all, it is worth noticing that for any function h £ L 1 (R), a > and r/Owe have 

\v\ a h\\ L i. (46) 



1 i ck 1 ? f V \ i 1 iqim 1 a j 

\v\ —fi y — J civ — \t\ 



Let us compute now the second moment of y^+i 

r- + oo 



/ 



v 2 tpf +1 (v)dv 



r 
At 



At 



- (<P?,p*<pL) (7) d ' 



v + 



+ (1 - At) 



« ^- ) dv 



dr 



(At\\v 2 ip* p * ip* q \\ L i + (1 - At)\\v 2 <p? |Ui) ^ 



T3T+ 1 

+°° r 2 dr 



7- At 



For AT > 221 we get 



dr 



At rit- 1 1 + At(p 2 

and so we are left with \\v 2 (pf p * tpj q \\ L i . We have 



1) 



2 1 N ( v-w\ 1 N ( w . 
K 2 p V P / 9 V " 



((« — w) 2 + w 2 + 2(w — iu)to) i</? 



V - W \ 1 N (W . 

p J q \q. 



Thanks to lf4"5]l. we get 



and we end up with 



11 2 JV JV n , 2 . 2mi 2 JVn 

I! 1 ' Vi,p * = (P +9 )||w ifij \\ L i 



11 2 at I, (At(p +9 -1) + 1),, 2 JV,, I, 2 JV,, 



I + At(p 2 + q 2 - 1) 
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and so by a recursive procedure 

[ v 2 <p? +1 (v)dv = l, j = 0,...,N-l. (47) 
As for J R \v \ 2+s ipf +1 {v) dv, we proceed in the same way: 

/ M 2+, Vj+i(«)du 

Jm. 



r 

At A 



= Jl (At|||H 2+i ^ * + (i - At)||M 2+i |Ux) ^ 

Now, for iV > ( 2 +^) T we get 



dr 

IF 1 

r 2+a d7 



dr 



At A ri.- 1 -' l-Al( 5 + 2 ) 1- At^±2(l-p 2 -g 2 )' 
Let us estimate |||i>| 2+lS </3.^ r p * y^gllii. We will denote 

d f = / \v\ 2+ "<pf(v)dv. 

Since < S < 1, we can write 

|«| 2+l5 _ w 2 |«;|* = (u — to -)- w) 2 \v — w + w| 5 < ((l> — ™) 2 + ™ 2 + 2(u — Ml) — to] 15 + |w|*J 

= \v — w\ 2+s + w 2 \v — w\ s + 2(w — w)\v — w\ S w + (v — to) 2 ]™] 15 + to] 2 ^ + 2(v — w)\w\ s w 



so, thanks to (|45[) and (|47|) . we get 



d,- +1 < 



1- At^(l-p 2 -g 



[dj + At(p +q T - l)dj + At(p q +q p )\\\v\ (fj \\ L i ) 



By Holder inequality and the conservation of the mass, we obtain 

\\\v\ 5 i pf\\ Ll <\\\v\ 2 <pf\\l 1 \\ V ,f\ff=l 
and so we get the recursive estimate 



1_ Ati±2(l-p2_ 



,N ^ ,N , A ./ 2+i . 2+d> i\jiV . a ,/ 4 2 4 2 

d j+i < ; A ,4+2/-, o — ( d i + At(p +9 " + At(p q +q p 



Remembering that At = we would like to neglect the low order terms. We have 
1- At^l-p 2 ^ 2 ) \ j + iP + q ~ ] 3 + (P q +QP 



2 

l + At^(l-p 2 -Q 2 j / 2+ , 2+ , 



(df + At(p 2+i + q 2+S - l)df + At(pV + q 5 p 2 ) 



l-(At) 2 (^(l-p 2 -g 2 )) i 
_ df + At(p 2+ * + g 2+S - 1 + *±3 (1 - p 2 - g 2 ))rif + At(pV + gV) + (Af) 2 JTdf + (AtfH 

l_(At) 2 (*± 2 (l-p 2 -g 2 )) 2 

where i/, if are positive constants, depending only on p, q and 5. Denoting 

<Sp, 9 (<5) = p 2+a + g 2+A - 1 + ^y^(! - p 2 - <7 2 ), (<*) = pV + <?V > 

we got 



< 



df + AtS p , q (S)df + AtB p , q (S) + [At) 2 Hdf + (At) K 



l-(At) 2 (^2(l-p 2 -g 2 ))' 
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where, by assumption, we have S Ptq (S) < 0. Moreover 

1 



l_ (At) 2(i±2 (1 _ p 2_ g2) ) 2 

so we get, for N £ N large enough, j = 0, . . . , N — 1 



l + o(Ai), At^O, 



<+i < < - ^At\S p , q (S)\d? + 2AtB p , q {5). 



f 



This recursive relation implies 

>\ 2+s V? +1 (v)dv<C s , j = 0,...,N-l, 
for any N £ N large enough. 

□ 

Proof of Proposition 1101 We divide it into several steps. 

I STEP: Existence of the limit <?*(£,£) of a subsequence. Thanks to inequalities (|26|) 
and to the definition of g N we have therefore 

sup \g N {Z,t)\<l, sup \d ( g N (Z,t)\<C, sup \d!g N (£,t)\ < 1, (48) 

[O.TJxR 1 1 [0,T]xR 1 [0,T]xR' 

where C is the same constant as in (I26p and iV is large enough. Moreover since g N (£,t) 
satisfies 

drg N (t,t) = ^^ w (o + ^-i(pe)^-i(90-^-i(o (49) 

then for any compact set i(cR there exists a constant C > such that 

I I 2T 1 

sup \drg N (t,t)\<C, N>—. (50) 

[0,T]xA- 1 1 r 

For any compact K C R the function g N ((,, t) belongs to C([0, T] x 7f) and thanks to properties 
(|48[) and (|50[) the sequence is equibounded and equicontinuous. Therefore by Ascoli-Arzela 
theorem and by taking the diagonal, there exists a subsequence {g Nl {(,,t)} which converges 
uniformly on [0, T] x K for any compact A' C R. Let us call <?*(£, i) the limit function. 
Since g N '(^,0) = ffo(C) then g*(£,0) = go{£.)- Moreover for any t £ [0, T] the function 
g*{£,t) £ C J (R): indeed by Ascoli-Arzela theorem and the diagonal argument applied now to 
the sequence {d^g Nl (£, t)}, for any t £ [0, T] we get a subsequence that converges uniformly 
to d^g*(^,t) on any compact set K C R. Since the limit function is d^g*{S,,t), the convergence 
holds for the original sequence {d^g Nl (£, i)} and it is not necessary to pass to a subsequence. 

In order to get a uniform convergence in both frequency and time, we remark that by (|48j) 
we have that supr T j xR l^lfir^^, i)| < 1 and thanks to (|49[) and (J26J) we control dtd^g N {(,,t) 
therefore d^g N ((,,t) is Lipschitz continuous on [0,T] x K and uniformely bounded. Again 
by Ascoli-Arzela we prove that {d^g Nl (£, t)} converges uniformly to d^g*(^,t) on [0, T] x K 
for any compact set K C R. From the uniform convergence of g N, (£„t) to p*(£,t) and of 
9ig N '(^,t) to d^g*(£,t) we get that both ff*(£,t) and <9 5 sf*(£,i) belong toC([0,T] x if). 



II STEP: g*(£,,t) is a solution of equation (|20[) . By a direct computation we obtain: 

Offf"' & *) = fleff^ K, *) + 9 Nl (p€, t)g Nl («£, t) - <? N < (£, t) + i?iv, «, t) (51) 

where 

(1 - a) [a (qO - (<fiZ, 0*) " + 



+ 



(52) 
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Since {Rn 1 (£, t)} converges to zero for any t £ [0, T] uniformely on any compact set fcl, 
the whole right hand side of equation (|5ip converges; let us call L(f , i) the limit function, so 
that 

= -eaes*«,*)+i7*(pe,«)s*(«e,*)-s*tt,*) (53) 

and {d^g Nl (£, t) } converges to L(£, i) for any t £ [0, T] uniformely on any compact set K C R. 
Thanks to property (|50p and to Lebesgue's dominated convergence theorem we remark that 
{dr9 Nl (£,t)} converges to L(£,t) in V'(]0,T[xK°). Since {g Nl (£,t)} converges to g*(£,t) 
uniformly on [0,T] x K and so in V'(]0,T[xK°) we have that {d t g Nl (f, t)} converges in 
distributions to &<;*(£, i). For the uniqueness of the limit and the fact that in distributions 
g Nl (£, t) = dtg Nl (£,t), we obtain dtg*(£,,t) = L(£,t) in the sense of distributions. Finally 
since g*(£,,t) and d^g*(^,t) belongs to C([0,T] x A"), for any compact K the right hand side 
of (|53[) is continuous in both variables therefore the same is true for dtg*(£,t). This implies 
that g*(£,t) belongs to C'([0,T] x R). 

Ill STEP: we show that g*(£,t) = fl(£,t). Defining F*(£,t) = g*(£y/E(t),t) and /(£,t) = 
g(£y/E(t), t) we obtain two solutions of equation (fT5)> . Both F* t) and /(^, are C'([0, T] x 
R) and bounded by 1. Denoting h(£,t) = f(£,t) — F*(£,t) by an easy computation we are 
led to 

\d t h + h\ < 2]|/i(i)||oo. 

By Gronwall Lemma we have 

l|ft(t)ll« <e*||ft(0)ll». 
and since /(^,0) = F*(£, 0) this implies the desired equality g*(£,,t) = g((,,t). 

□ 
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